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Abstract: A new method is presented for measuring the spin of selectrons and 
smuons at the Large Hadron Collider (LHC), using an angular variable which is sen- 
sitive to the polar angle in direct slepton pair production. This variable, cos 6^, is 
invariant under boosts along the beam axis, so it can be used at the LHC despite 
the fact that the longitudinal boost of the centre-of-mass frame cannot be deter- 
mined. Monte Carlo simulations demonstrate that, using this method, the LHC can 
distinguish between the supersymmetric production angular distribution and phase 
space, or between supersymmetry and the production angular distribution of uni- 
versal extra dimensions. An integrated luminosity of about 100 to 300 fb _1 provides 
sufficient statistics to measure the slepton spin for points which had left-handed 
slepton masses in the range 202 to 338 GeV, and right-handed sleptons in the range 
143 to 252 GeV. Good sensitivity was found in the 'bulk' and 'stau co-annihilation' 
regions of the cMSSM favoured by cosmological relic density measurements. Various 
systematic uncertainties are investigated, and some methods for reducing them are 
discussed. 
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1. Introduction 

One of the main goals of the LHC is to search for new physics, beyond the Standard 
Model (SM). In most cases, if new particles exist at the TeV-scale, then the LHC 
experiments will be able to discover their existence. But discovery is only the first 
step - to distinguish between models, the properties of the new particles also need 
to be experimentally determined. 

Probably the leading contender for new physics at the TeV scale is supersymme- 
try (SUSY). The minimal supersymmetric extension to the SM predicts that each SM 
particle should have a partner with spin differing by |. If R-parity is conserved, then 
the supersymmetric partners must be pair-produced and the lightest supersymmetric 
particle (LSP), typically the %\i i s stable. The typical LHC signature is therefore 
a pair of cascade decays which produce jets and leptons, as well as missing energy 
from the (assumed invisible) lightest supersymmetric particle. 

There are, however many other competing models for new physics at the TeV 
scale besides supersymmetry. Indeed it has recently been appreciated that models 
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with Universal Extra Dimensions (UED) and Kaluza-Klein (KK) parity [1] could 
have very similar collider phenomenology to supersymmetric model [2]. 

The minimal version of UED predicts that for each SM particle there should be 
a tower of Kaluza-Klein (KK) excitations. KK parity means that particles with odd 
KK-number, such as the first excited state of any SM particle can only be produced 
in pairs. It also ensures that the lightest KK particle must be stable, in the same way 
as R-parity does for supersymmetry. Distinguishing between SUSY and UED could 
therefore be a difficult problem, since both models predict TeV-scale pair-produced 
particles which decay through cascades with Standard Model couplings, with the 
eventual production of a pair of invisible daughters (LSP or LKP). 

While other measurements might be indicative [3,4], the property which will a 
give conclusive answer as to whether we are observing SUSY or UED is the spin of 
the excited particles. 

Recently some progress has been made towards spin-determination of supersym- 
metric particles at the LHC. The method, suggested in [5] and investigated in [3,6,7], 
involved measurement of the lepton charge asymmetry in £q invariant mass distribu- 
tions in the cascade decay, 

QL^xUL^liF QL- (1-1) 

That measurement was shown to have sensitivity to the spin of the x\- While it was 
comforting to see that the LHC can have sensitivity to sparticle spins, the caveat 
is that in some parts of parameter space, the decay chain eq. is kinematically 
forbidden or has a small branching ratio. This makes it important to investigate 
other channels and other particles for which the LHC experiments could measure 
spin. 

In this paper we present a new method for measuring slepton spin at the LHC. 
The paper is organised as follows. In section [| we introduce an angular variable 
cos 6^, and show that it is sensitive to the production polar angle in slepton pair pro- 
duction. Our supersymmetric test points, Monte Carlo event generator and detector 



simulation are described in section |3|. In section |3.1| we identify an event selec- 
tion and demonstrate that it can cleanly isolate the signal process. Results showing 
the experimentally-measurable angular distributions and luminosity requirements are 
shown in section [3.2| . In section || we discuss the main systematic uncertainties and 
some methods for reducing them. Our conclusions are presented in section ^. 

2. Angular distributions, and cos 6^ 

In this paper we investigate the supersymmetric process, 

qq - Z*h - - X\t X?r , (2.1) 
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where throughout this paper t is understood to mean electron or muon only. Since 
sleptons are scalars, the angluar distribution for Drell-Yan slepton pair production 
is 

' da \ ., 

ocl-cos 2 r (2.2) 



d cos 9* 



SUSY 



where 9* is the angle between the incoming quark in one of the protons and the pro- 
duced slepton. Slepton pair production via gauge boson fusion [8] is not considered 
here, but it would become important for sleptons with masses greater than about 
300 — 400 GeV. For comparison we use a pure phase space distribution, 



da 



d cos 9* 



oc constant 



(2.3) 



PS 



The phase space distribution does not correspond to any physical model, but does 
provide a convenient benchmark against which to compare the SUSY distribution. 
We also compare against the UED equiv- 



alent of eq. 2.1 



qq 



z°h - *t ^ - 7i £ + n r . (2.4) 

which has the characteristic distribution 
for spin-- KK leptons: 



da 



d cos 9* 



OC l+i 



VED 



cos 2 9* , 
(2.5) 
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where and are the energy and mass 
respectively of the KK leptons in the center- 
of-mass frame. The three different pro- Figure 1: 
duction angular distributions are shown grapfi lons ' icose 
ically in fig [I] spin-^KK leptons UED and pure phase 

The different angular distributions pro- s P ace ( PS )- The mass s P ectrum for the 



Production angular distribu- 
-, for scalar sleptons (SUSY), 



UED distribution is that of SUSY point S5 
(see section 



vide a mechanism for determining the heavy 
particle spin. Excited leptons (selptons or 
KK-leptons) which are produced significantly above threshold will have decays which 
are boosted in the lab frame. This means that a pair of leptons from slepton decays 
(eq. [2.2|) should be on average less widely separated in polar angle than the pair from 
phase space (eq. |2T3| ) or KK-lepton pair production (eq. |2~5| ). 

It has already been suggested [9, 10] that the final state lepton angular distri- 
butions could be used at a future high-energy e + e~ linear collider to distinguish 
between UED and SUSY models. With a proton-proton collider such as the LHC, 
it is not possible to measure the lepton angluar distributions in the parton-parton 
center-of-mass frame - the initial ^-momenta of the incoming partons are not known, 
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Figure 2: 2-dimensional plots showing the correlation between our dilepton angular vari- 
able, cosfy*, (y-axes) and the cosine of the production angle of the parent sleptons (a) or 
KK-leptons (b) in the center of mass frame (x-axes). Darker regions correspond to larger 
numbers of events, the normalisation being arbitrary. The mass spectrum is that of SUSY 
point S5. 

and because invisible particles are produced, the center of mass frame of the parton- 
parton interaction cannot be recovered from the final state. 

To make a spin measurement at a hadron collider, we propose a variable which 
is a function only of the pseudorapidity difference between the final state leptons, 
Ai]i+i-. The advantage of differences in rapidity is that they are independant of the 
longitudinal boost. The leptons are highly relativistic, so we can use their pseudo- 
rapidities as a very good approximation to their true rapidities. By using a function 
only of Ar]£+£-, we no longer have to face the problem of determining the center- 
of-mass frame along the beam direction. The inter-lepton pseudorapdity difference, 
Ar]£+£~, is also sensitive to the slepton production angle. The reasons are the same 
reasons as for the lepton angular distributions - the leptons 'inherit' some knowledge 
of the rapidity of their slepton or KK-lepton parents. Lepton pairs from slepton pair 
decay will therefore be on average less separated in pseudorapidty than those coming 
from particles produced according to the corresponding phase-space or Kaluza Klein 
production angular distributions. 

To allow a more direct comparison with the production distributions, rather than 
using Ai]£+£- directly, we propose the angular variable 

cos^ = cos ^2tan _1 exp(Ar/£+£-/2)^ = tanh(A?7£+£-/2) . (2.6) 

This variable, like At)i+£-, has the benefit of being longitudinally boost invariant, 
but also has a simpler geometrical interpretation: cosfy* is the cosine of the polar 
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angle between each lepton and the beam axis in the longitudinally boosted frame in 
which the (pseudo-)rapidities of the leptons are equal and opposite. 

As can be seen from fig. [^, cos 6?^ is indeed well correlated with the slepton or 
KK-lepton production angle, cos 9*. The experimental observable cos 9^ is on average 
smaller for SUSY (fig. |]a) than for UED (fig. |2|b), meaning that cos 9* t can be usefully 
employed as a spin-sensitive discriminant in slepton/KK-lepton pair production at 
hadron colliders. 

3. Monte Carlo simulation of test points 

We chose as a primary test point the LHC point 5 [11], (S5), with supersymmetric 
mass spectrum and branching ratios calculated with ISAJET-7.64. Further simu- 
lations were performed with the Snowmass points SPSla, SPSlb, SPS2, SPS3, SPS4, 
SPS5, and SPS6 [12] with spectra and branching ratios calculated using ISAJET-7 . 58, 
following the Snowmass standard. 

Some important masses as well as the signal cross-section are given in table [TJ 
It is immediately clear that SPS2, which is in the 'focus point' 1 region, and has very 
heavy sleptons, will not provide a measurable signal at the LHC. Of the other points, 
we note that SPS4 has slepton masses near the edge of our expected sensitivity for 
direct production [14]. The non- universal point SPS6 has a smaller-than- usual mass 
difference between the sleptons and the LSP, and so will produce rather soft leptons 
in its decays. 

The Monte Carlo event generator used was HERWIG-6 . 507 [15-17], with the par- 
ton distributions of MRST [18] (average of central and higher gluon). Since a complete 

1 For a discussion of the various regions of cMSSM consistent with dark matter see, for 
example, [13]. 



Point 


X? 


en 




a(pp g+f) 


S5 


116 


153 


229 


65 fb 


SPSla 


96 


143 


202 


88 fb 


SPSlb 


160 


252 


338 


12 fb 


SPS2 


80 


1452 


1456 


0.003 fb 


SPS3 


161 


178 


287 


34 fb 


SPS4 


119 


417 


448 


2.5 fb 


SPS5 


120 


191 


256 


48 fb 


SPS6 


189 


237 


265 


22 fb 



Table 1: Masses of selected particles (in GeV) for the model points investigated, and the 
sum of the leading order cross-sections for direct slepton (e + e~ and pair production. 
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UED Monte Carlo event generator programme was not available, HERWIG was modi- 
fied to generate slepton pairs with the other production angular distributions, eq. p73| 
and p.5| . The slepton angular distribution (eq. |2.2|) is of course already implemented 



in HERWIG [16]. Angular distributions for UED processes other than for t\l\ pro- 
duction in eq. p.4| were not modified. In general, this should not significantly affect 



the results presented, since cuts were later applied (section |3.1| ) which reduced the 
residual supersymmetric background to a small fraction of the signal. 

Unlike the case of slepton decay, the extent to which the spin structure of the 
^1 ~~ > vertex is important will depend on the details (such as the KK Weinberg 
angle and left-right i\ mixings) of the UED model. We expect that others with an 
interest in particular UED models will investigate in more detail the performance of 
this method for those models. For this paper, we merely note that the KK-leptons 
are only weakly polarised (as can be seen from the similarity between the UED and 
PS distributions in fig. |I|), and so in our simple approximation to UED we ignore 
spin correlations in the t\ decays. 

Samples of inclusive supersymmetric particle production were generated corre- 
sponding to integrated luminosities from 25 to 500 fb" 1 for each of the test points. 
Also, 500 fb _1 samples were generated for the major backgrounds, iy ± iy =F , Z°Z°, 
W ± Z°, and ti. To understand the 'theoretical' distributions, very large luminosity 
(10 ab _1 ), samples were generated of the signal process for each of the three produc- 
tion angular distributions (eqs. p72| , j2~3[ and |2.5|) , and for each test point. In total 



about 8 x 10 8 events were generated. 

All of the events generated were passed through the ATLAS fast detector sim- 
ulation programme ATLFAST-2.50 [19], which contains a parameterisation of the 
ATLAS detector response for leptons, jets and missing energy. Jets were defined 
by having 10 GeV of transverse energy in a cone of radius 0.4. The pseudorapidity 
coverage for jets was \rjj\ < 5, while the lepton {t e e,/i) pseudorapidity coverage 
was \r)g\ < 2.5. We assume a 6-jet vertex tagging efficiency of 60%, a c-jet vertex 
tagging rate of 10% and light-quark jet mistagging rate of 1% [11]. 

3.1 Event selection 

The final state of interest consists of two opposite-sign, same family (OSSF) leptons, 
either electrons or muons, together with two invisible particles. Thus the initial 
selection is a requirement that exactly two OSSF leptons are detected, and that there 
must be missing transverse momentum, fa > 100 GeV. The lepton with the larger 
transverse momentum must satisfy pr(^ ) > 40 GeV, while that with the smaller 
Pt must have pr(^) > 30 GeV. This lepton momentum selection also ensures that 
these events satisfy the ATLAS trigger requrements. 

The main Standard Model backgrounds are expected to come from the final 
state i^i'vv, from H /± VT =F or Z°Z° production. Since W^Z and ti might also be 
significant, they were also generated. 
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The backgrounds is largely removed by excluding events with dilepton 

invariant mass mu. < 150 GeV. 

The W /± W' r=F pair background can be substantially reduced by making a cut on 
a different kinematic variable, known as m^. This variable, described in [20,21], 
is similar to the transverse mass, tut, but is useful in events in which a pair of 
same-mass particles decays semi-invisibly. It is defined by: 



mT2 2 (mu) = um) 



max {m^p^ , $r ; m v ) , m|(p^ 2> , <f$\ m ")} ■ 



m 2 T (p e T , p v T - m u ) = m]+ + m 2 u + 2[E l T E v T - p e T ■ p v T ) , (3.2) 



where, 



Et = V(Pt) 2 + mj and E U T = ^(p^) 2 + m 2 . (3.3) 

and p^. and p^>, indicate the lepton and neutrino transverse momenta respectively 2 , 
mn has the property by construction that for events in which the lepton pair origi- 
nates from on- mass-shell W ± W T pairs, 

m T 2(PT,PT 'i m v) < m w ■ (3.4) 

We therefore exclude from our analysis any event for which 771t2(Pt )Pt > 0) < 
100 GeV, removing most of the W /± W /=F background. 

To reduce the background from ti production and from cascade decays of heavier 
KK or SUSY particles, the following additional selection was applied: 

• no jet with p T > 100 GeV 

• no jet with a vertex tag (no '6-jets') 

• transverse recoil, |pV + Pr(^ + ) + Pt(^~)| < 100 GeV, 

where the jet and vertex tagging algorithm parameters are as described previously. 

Another potential background comes from single W ± production in association 
with a jet or photon which fakes a lepton in the detector. The evaluation of such 
backgrounds would require a complete detector description, beyond the scope of this 
paper, and so they have not been included in this analyis. In fact we expect that 
although they have a much higher cross-section than gauge boson pair production, 
they should not pose a serious threat. The expected rate of jets faking electrons is 
rather small - of the order 10 -4 — 10~ 5 [11]. What is more, these types of event 
should also be efficiently removed by the cuts. Fake events where the is close to 
its mass shell and in addition either: 

2 The neutrino mass is not important in this analysis, but is included in eqs. 3.1 ■ 3.3] for com- 
pleteness. 
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• there is no significant missing transverse momentum (j^r) except that caused 
by the neutrino from the W ± decay or 

• any additional contribution to j/>t is collinear with the 'fake' lepton, as could 
be expected from a poor energy measurement of that fake, 

will fail the mn cut. The reason can be understood when one considers the quantity 
being minimized in eq. |3.1| with a particular pair of value for the 2-vectors $p and 

(2) 

$p which represent a possible transverse missing energy spitting. The interesting 
configuration is one in which $p is equal to the true neutrino transvers momentum, 
and q^ is parallel to the fake lepton. Under the conditions above, this configura- 
tion will satisfy the constraint + = The first transeverse mass, mS^ , 
containing the neutrino and lepton from the W ± decay, must be less than mw, 

(2) 

while the other, m T , containing the fake lepton, will be close to zero. Since such a 
configuration exists, it places an upper limit on of raw- 




Figure 3: Distribution of cos 9*i for the S5 SUSY sample after event selection. The dark 
cicles represents the SUSY slepton pair signal (eq. ^3). The open cicles indicate the SUSY 
background. The open squares show the different flavour (OSDF) estimate of the SUSY 



background, which is explained in section 4.4. In this and in all subsequent plots, the 



electron-pair and muon-pair distributions have been combined. 

The cuts are somewhat similar to those used for studies of slepton discovery 
potential, or mass measurements from Drell-Yan pair production at the LHC [14, 
22-26]. One notable difference is that, unlike most previous studies, we make no 
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Figure 4: The Standard Model backgrounds after the cuts have been applied, as a function 
of cos Of,. 



azimuthal angle (0) requirements on the leptons or ji r vector, since this seem to be 
unnecessary when an m T2 cut is used to reduce the W ± W T background. 

Fig. | shows that after the cuts have been applied, the SUSY background is 
very much smaller than the signal. Similarly fig. |] demonstrates that the Standard 
Model backgrounds do not overwhelm the SUSY signal. However the SM background 
is not insignificant, and so it will have to be accurately determined if a slepton 
spin measurement based on cos 6^ is to be achieved. The degree to which these 



backgrounds can be understood is investigated further in sections [O] and \L3 



3.2 Results 

In this section we disuss the cos 9* t distributions which would be obtained assuming 
that the supersymmetric and Standard Model backgrounds can be accurately deter- 
mined. We then go on to examine the extent to which these assumptions can be 
justified in section |[ 

In fig. [5] we plot the fraction of events as a function of cosfljlj for 200 fb -1 of 
integrated luminosity. As has already been shown in figs. [| and ||, for point S5 the 
statistical uncertainty on the SUSY signal and the SM background are comparable, 
while the SUSY background is small. It is clear that the (SUSY) "data" sample is 
much better matched to the slepton angular distribution than to either the phase- 
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Figure 5: The points show the cos 0T, distribution for the S5 signal sample 



after an integrated luminosity of 200 fb . The lines show the predictions for 



angular distributions according to supersymmetry (solid black line, eq. 2.2), phase space 
(dotted blue line, eq. ^1 ), and universal extra dimensions (dashed red line, eq. |2.5| ). The 
error bars on the data show the statistical uncerainty on: inner error bar: SUSY signal 
only; intermediate error bar: inclusive SUSY with the SUSY background subtracted; outer 
error: inclusive SUSY with both the SUSY and the SM backgrounds subtracted. The 
narrow shaded band around the SUSY expectation shows how it is modified when the 
sparticle masses are simultaneously changed for all sparticles by ±20 GeV, as described in 
section 4.4. Systematic uncertainties in the SUSY and SM background subtraction are not 
included here, but are discussed in sections 0~2| and [O. 



space one or the UED-like one. This means cos 6^ does indeed measure the spin of 
the sleptons for this point. 

In fig. | we present the statistical separation expected for our test points (S5 
and the Snowmass points) as a function of integrated luminosity. The significance 
indicated is shows the gaussian-equivalent significance of each of two tests: 

1. A test comparing the SUSY angular distribution (eq. [2~2] ) to the phase space 
one (eq. |2.3| ) - demonstrating that there is sensitivity to spin in the dynamics; 

and separately, 

2. A test comparing the SUSY angular distribution to the UED-like one (eq. |2~5| ) 
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Figure 6: The expected statistical significance for discriminating SUSY from UED (solid 
lines), or SUSY from phase space (dashed lines) for our SUSY test points, as a function of 
integrated luminsity. Statistical uncertainties are included for the signal, the SUSY back- 
ground from slepton-pair to gaugino (other than Xi) an d the Standard Model backgrounds 
in fig. ||. For each test point the the solid line (SUSY vs UED) always requires lower 
luminosity. 

- showing discrimination between two physically-motivated models. 

In both cases we have used a discriminant which is symmetrical in the hypotheses 
under test. The discriminant accounts for the fact that, because the normalisation of 
the distribution is fixed, the data in the bins are correlated. The events were counted 
in large statistics samples after the cuts had been applied and include the statistical 
uncertainty from the SUSY and SM background and the slepton to gaugino (other 
than Xi) background subtraction, but none of the other systematic uncertainties 
discussed in section f|. 

Starting with the worst case, SPS2 (from the cosmological 'focus point' region) 
has > 1 TeV mass sleptons, and such a small cross-section for direct slepton pro- 
duction that no spin measurement (or indeed any other direct slepton production 
measurement) is possible for this point. It is not included in fig. |6|. 

SPS6 (also not on fig. |6[) also presents a difficult experimental case, but for a 
quite different reason. This point has non-universal gaugino masses, with the x° 
more massive than would be the case with universal gauginos. Despite its large 
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(d) 



Figure 7: As for fig. || but for Snowmass points (a) SPSla, (b) SPSlb (c) SPS3 and (d) 
SPS5. The integrated luminosity simulated in these plots is 200 fb -1 for SPSlb and SPS 3, 
and 300 ftr 1 for SPS la and SPS 5. 



cross-section, a spin determination is very difficult because the sleptons are only 
about 50, and 75 GeV heavier than the LSP, and the low-pr leptons they produce 
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end up buried beneath the W /± H^ =F background. 

SPS4 has tan/3 = 50, rather heavy e and jl sleptons (table [I]) but relatively 
light xf an d X2 a ^ 218 GeV. The effect of this is that the signal process eq. |2TT] has 
to compete with cascade decays through the xf or X2 anc ^ so the super symmetric 
background is rather large at this point. For this point, either very high statistics 
(such as could be achieved with an LHC luminosity upgrade [27]) or perhaps a very 
careful background reduction analysis, might allow the spin determination to be 
made in this channel. 

The remaining five points, S5, SPSla, SPSlb, SPS3, and SPS5 all allow spin 
determination with 100 — 300 ftT 1 of integrated luminosity. Reference distributions, 
similar to fig. [5], but for SPS points, are shown in fig. [?| 

S5 and SPSlb which have masses 'typical' of the SUSY 'bulk' region require 
integrated luminosity of about 100 and 300 fb _1 respectively. The point SPS5, which 
is distinguished by a very light stop squark, 220 GeV, behaves in this analysis like any 
other 'typical' spectrum. SPSla requires a slightly higher than expected integrated 
luminosity (about 300 fb _1 ) despite it having the largest production cross-section 
(table p. This is because its lighter sleptons produce relatively soft leptons which 
more frequently fail the cuts, and so the cross-section after selection is actually 
smaller than for many of the other points. 

Point SPS3 is in the cosmo logical 'stau co- annihilation' region, in which the 
lighter stau is nearly degenerate with the Xv One might worry that the leptons it 
would produce might be too soft for this analysis. While this is true for the right- 
handed sleptons, which are only 17 GeV heavier than the LSP, a good signal is found 
from the heavier left-handed sleptons, and it turns out that this point is in fact one 
of the easier ones with which to make this spin measurement. 

In one sense fig. |B| is optimistic, in that it does not correctly take account of all 
the experimental systematic uncertainties explored in section |] - only the statistical 
effect of the SUSY signal, and the SUSY di-slepton to gaugino and Standard Model 
background subtractions are included. In another sense this plot is conservative, as 
no tuning of the cuts has been done between the different points. With the real data, 
the cuts could be tuned according to the masses of the particles already measured. 
Such tunings have not been attempted here but have previously been shown to pro- 
vide somewhat improved results for making mass measurements at different slepton 
masses [14,24]. 

Overall, we can see that, unlike the previous lepton charge asymmetry mea- 
surement [5], our new method appears to be rather general. There are some difficult 
cases, when the sleptons are heavy or the slepton-LSP mass difference is small. How- 
ever spin determination is reachable with statistics corresponding to a few years of 
LHC design luminosity for a variety of points coming from both the 'bulk' and the 
'stau co-annihilation' regions favoured by cosmology. 
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4. Systematic uncertainties 

It is important that neither the experimental data points nor the comparison distri- 
butions are biased in a way which cannot be either measured or calculated. Having 
shown in section |3] that the spin measurement is within statistical reach, we now 
address the question of how the various systematic uncertainties can be controlled. 

4.1 Lepton efficiency, acceptance and charge identification 




Figure 8: Reconstructed pseudorapidity (7/1,172) of the each of the two leptons in a large 
sample of signal events at SUSY point S5. The ATLFAST lepton acceptance limits at |t/| = 
2.5 are indicated by the dashed lines. 

Since the sensitive variable is based on At]£+£-, it is clear that the electron 
and muon reconstruction efficiency will have to be well known across the range of 
pseudorapidity. In ATLAS, the reconstruction efficiency is expected to be accu- 
rately determined from the experimental data, using the high statistics sample of 
Z° — > available at the LHC. One of the leptons is used as a tag and the 

efficiency is determined from the fraction of events in which the second lepton is 
identified. The same sample, and the same methods can be used for determining the 
charge mis-identification rate. This can be done for different values of rj, and so the 
reconstruction efficiency should be well under control for all values of 77, including 
around the barrel-endcap transition, and near the limits of the detector acceptance. 

Focusing our attention more closely on the edge of the detector acceptance, we 
note that it is partly due loss of events beyond acceptance that there is a down-turn 



14 



□ WW enhanced 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



cos 8n 



Figure 9: cos fy* distribution for a control sample in which the W^W^ pair contribution 
has been enhanced by reversing the cut on rriT2- 

in the cos Q* n distributions (figs. || and |7|) near unity. To investigate the sensitivity 
to the experimentally-forced rapidity selection we plot in fig. || the pseudorapidity 
distributions of the selected leptons. Only a small proportion of entries occur near 
around (|?7| = 2.5), suggesting that uncertainties in the acceptance should not cause 
a significant bias in the results. 

4.2 Standard Model background determination 

The normalisation for each Standard Model background can be determined indepen- 
dantly from the data, by measuring the cross-section in regions of phase space in 
which each particular background dominates. For example, if one reverses the 
cut, so that only events with m,T2 < 100 GeV are accepted, one obtains a sample 
of events rich in W ± W Zf boson pairs (fig. |9|). In the same spirit, it is not difficult 
to define cuts which preferentially select Z° boson pairs, W ± Z° events or tt pairs. 
Some care obviously needs to be taken to ensure that a bias is not introduced when 
extrapolating from one kinematical region to another. The experimental techniques 
required for such an extrapolation are not investigated here, but will be the subject of 
a future study [28]. Previous work on SM background subtraction using the flavour 



subtraction methods discussed in section fO] has already shown good promise [14] . 
4.3 SUSY backgrounds 

Although the SUSY backgrounds are generally small compared to the signal (fig. |3|), 
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they do need to be understood so that they can either be confirmed to be too small 
to worry about, or, more realistically, so that they can be accurately determined and 
subtracted. 

Frequently-used methods [7,24] for measuring the background rate from the data 
are based on constructing the equivalent distribution as for the signal, but rather than 
requiring an opposite sign, same family (OSSF) dilepton pair, asking instead for a 
same sign, same family (SSSF e.g. = e + ,e + ) pair; or an opposite sign, different 
family (OSDF, e.g. /i + e~) pair; or a same sign, different family (SSDF e.g. e + n + ) 
pair. Combinations of SSSF, OSDF, and SSDF selections can be employed along 
with particular background-enhanced samples to estimate the various backgrounds. 

For high-scale processes at the LHC a flavour subtraction rather than a charge 
subtraction method is generally most powerful. Lepton universality is only broken by 
the difference between the \x and the e Yukawa couplings, so distributions involving 
hard, isolated leptons (as opposed to those coming from near-threshold or Dalitz 
decays) can be generally expected to be rather universal between generations. By 
contrast at high scales we are sampling the proton's parton distribution functions in 
the valence region, so the initial state is not charge symmetric, and so neither are 
the final states, reducing the effectiveness of a charge-subtraction method. 

For slepton pair production most of the SUSY backgrounds can be expected to 
have similar cross-sections in the OSSF (signal) and the OSDF (control). The good 
similarity between the SUSY background and an OSDF selection for our S5 point 
can be seen in fig. |3|. It has been shown [14] that a simple OSSF-OSDF subtraction 
(without using control samples) can provide a rather clean background subtraction 
for both the SUSY background and most of the Standard Model backgrounds. 

Returning for a moment to the SM backgrounds discussed in section fO|, the 
size of the W ± W Zf , W ± Z°, ti, and + 'fake lepton' can also be estimated with a 
flavour subtraction. It is likely that the best results for both the SUSY and the SM 
background determinations will come from a mixture of techniques, both of OSDF 
subtraction and control samples. 



4.4 SUSY spectrum 

Another important issue to address is to what extent our lack of knowledge of the 
masses of the SUSY particles might lead to uncertainties in the comparison ('the- 
oretical') distributions. It is safe to assume that because of the higher integrated 
luminosity required, any spin measurement will be made after some knowledge of 
the sparticles masses has been obtained. It is therefore sufficient to investigate the 
effect on the comparison distributions caused by varying the sparticle masses within 
the limits which are expected after 100 fb _1 or more has been collected. 

Studies 3 suggest that with around 100 fb _1 of integrated luminosity the LHC 

3 Examples include [11,25,29-31]. 
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experiments can strongly constrain the difference in masses between sparticles, but 
that the overall mass scale might only be known to a precision of about 10%. To 
determine the effect the sparticle spectrum has on the SUSY comparison distribution, 
the masses of all of the sparticles were simultaneously raised and lowered by 20 GeV. 
The change in the S5 comparison distribution, indicted by the shaded band in fig. |5] 
was small compared to the statistical uncertainty, so we can conclude that, for this 
point, lack of knowledge of the SUSY mass scale will not be a significant source of 
systematic uncertainty. A similar test was not performed for the other points, but it 
can be observed (fig. ^) that the shape of the cos 6^ distributions is rather universal, 
despite the differences in the spectra. 

4.5 Migrations at cut boundaries 




100 200 50 100 100 



p T (j (1) ) Pt^ij) recoil 




Figure 10: Distributions showing the number of events close to each of the various cuts, 
for signal (thick line) and backgrounds (thin lines, not to scale). The variables are the 
same selection as those described in section |3.1[ Note that the distributions have been 
made after the cuts have been applied. The position of the cut is indicated by the vertical 
dashed line. 



If the distribution of interest has a strong dependance on a poorly-known selec- 
tion variable, this can generate a systematic uncertainty. Migrations might occur, 
for example, from uncertainty in the hadronic energy scale, or from larger jet- or 
missing-energy resolution due to pile-up (minimum bias) events. In fig. [10] the num- 
ber of signal and SM background events (not to scale) are plotted as a funtion of 
the selection variables. The leading jet p?, leading lepton p?, and recoil have very 
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few events near their corresponding cut, so migrations due to uncertainties in these 
variables should not significantly affect the final distributions. 

The leading jet pt and recoil cuts are about a factor of two larger than the value 
(« 40 GeV) at which pile-up events are expected to cause fake vetos at the percent 
level [11]. The number of signal events lost by jets coming from pile-up events should 
therefore be negligible in this analysis. 

The three variables on the bottom line of fig. [H], fa, ttlt2, and mu, do have a 
significant fraction of events close to their respective cuts, so a more careful study 
of the resolution in these variables would be appropriate. For this analysis one is 
particularly interested in finding out to what extent migrations might change the 
shape of the cos 9f t distribution. 

While a detailed determination of the systematic uncertainties requires a more 
complete detector simulation, and is beyond the scope of this paper, we suggest that 
we might already have some optimism about the results. A large proportion of the 
visible transverse energy in the signal events is carried by the final state leptons, 
which should be well-measured in the electromagnetic calorimeter (electrons) or in 
the tracking system (muons). This means that the missing energy should also be 
relatively well measured, and since all three of the variables m«, fa, and m^, depend 
only on pt £ , and fa, they ought to be well under control. 

4.6 Migrations between bins 

Our bin size of A cos 9*i = 0.125 corresponds to pseudorapidity differences of Arji+i- > 
0.25. By contrast the r\ segmentation in the first sampling of the ATLAS liquid argon 
calorimeter in the \rj\ < 2.5 region is in the range 0.003 to 0.006, depending on rapid- 
ity. This fine r\ segmentation will provide a much better resolution than is required 
for our spin measurement. Similarly, for the muon reconstruction, the rapidity mea- 
surement from the combined inner tracker and muon chambers will be much more 
precise than is required here - inter-bin migrations should not cause problems. 

4.7 Summary of systemic uncertainties 

In this section we have explored the various experimental systematic uncertainties 
which could potentially reduce the power of the cos# z * z distribution in making a 
spin measurement. We have found that, while there are some areas - such as in 
SM background estimation from data, and in understanding the effect of the missing 
energy resolution - which warrant further attention, there is no reason to believe there 
are any serious problems - no 'show-stoppers' which could invalidate the method. 

5. Conclusions 

We have described a new method for slepton spin determination at hadron colliders. 
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It is based on a measurement of the di-slepton production anguar distribution, using 
a longitudinally boost-invariant variable we have called cosfy*. 

The method is rather general, and is able to measure the spin for variety of 
supersymmetric model points. 

We have examined LHC test point 5 (S5) , and seven of the Snowmass test points 
(SPSla— >SPS6). Generally the results were very encouraging, though some of the 
test points presented particular problems. The method was unable to make a spin 
measurement at SPS2 (focus point region) because the production cross-section was 
much too small, nor at SPS6 (non-universal point) because the difference in mass 
between the sleptons and the lightest neutralino is small. SPS4 is a difficult case, 
with relatively heavy sleptons, and a large SUSY backround. 

Each of the other five points tested - LHC point 5, SPSla, SPSlb, SPS3, and SPS 
5 - which include cosmological 'bulk' points, a 'stau co-annihilation' point, and a 
point with a light t\ allowed slepton spin determination with integrated luminosity 
in the range 100 to 300 fb~\ 

There do not appear to be any systematic uncertainties which might be seriously 
detrimental to the measurement, but we indicate that more work does need to be 
done in this area. 
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